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Abstract 

Supermassive black holes (BHs) appear to be ubiquitous at the center of all 
galaxies which have been observed at high enough sensitivities and resolution 
with the Hubble Space Telescope. Their masses are found to be tightly linked 
with the masses and velocity dispersions of their host galaxies. On the other 
hand, BHs are widely held to constitute the central engines of quasars and 
active galactic nuclei (AGN) in general. It is however still unclear how BHs 
have grown, and whether they have co-evolved with their hosts. In this 
Review I discuss how, in ways independent of specific models, constraints on 
the growth history of BHs and their host galaxies have been set by matching 
the statistics of local BHs to the emissivity, number density, and clustering 
properties of AGNs at different cosmological epochs. I also present some new 
results obtained through a novel numerical code which evolves the BH mass 
function and clustering adopting broad distributions of Eddington ratios. I 
finally review BH evolution in a wider cosmological context, connecting BH 
growth to galaxy evolution. 
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1. Introduction 

Supermassive black holes (BHs) are common at the center of most, if 
not all, galaxies (e.g., Magorrian et al., 1998; Ferrarese and Merritt, 2000; 
Gebhardt et al., 2000). As it will be discussed below, their mass is tightly 
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linked with several properties of their host galaxies. However, it still con- 
stitutes a major challenge understanding how these BHs have grown, and if 
they followed the mass assembly history of their host galaxies. 

The discovery of such strong correlations lent further sup port to the sem- 
inal idea, already p roposed more than forty years back (e.g., ISalpeterl . Il964 ; 



Lynden-Belll . 119691 ). that quasars and, more generally, active galactic nuclei 



(AGNs) are powered by mass accretion onto central BHs, which should now 
lie at the center of galaxies as relics of a past luminous phase. The only 
process in fact able to generate the high and highly concentrated powers 
observed in quasars is the release of gravitational energy from an infalling 
body of mass m onto an extremely compact object of mass M. This process 
can then generate energies up to a maximum ol E = GMm/5Rg = 0.1 mc^, 
with the Schwarzschild radius Rg = 2GM/c^ ~ 3 x 10^^(M/10^ M©) cm. 
Gravitational processes can therefore release significant fractions of the rest- 
mass energy, orders of magnitude more efficiently than the energy release 
from stellar nuclear reactions. For example, a quasar with luminosity of 
L = lO^^ergs"^ can be easily generated by a massive object of 10^ Mq con- 
fined within a region of size < 10^^ cm. 

The detection of any massive dark object at the center of a currently 
inactive galaxy is performed by analyzing the motion of gravitating bodies 
embedded in the strong BH gravitational field. In any given galaxy, a BH 
detection can be roughly labeled in terms of r^, the radius of the sphere 
of infiuence within which the BH dominates the gravitational potential of 
its host galaxy. However, it is easily seen that for_a 10*^ M©, at a distance 
of 20 Mpc, vg is only 0.11 arcsec (e.g., iBarthl . |200J), which proves that 
reliable direct dynamical BH mass measurements can only be safely obtained 
in nearby galaxies. In particular, the presence of BHs in inactive galaxies 
is te sted through the velocit y dispersion of stars, water masers, and gas 
(see 



Ferrarese and FordI l2005l ). The possibly most exciting example is the 



detection of a BH at the center of our own Galaxy, which also represents one 
of the tightest constraints on the existence of a BH. Monitoring has been 
carried out for more than 10 years on 40 stars orbiting around the center of 
the Galaxy. Measurements in the infrared with adaptive optics imaging of the 
10m class telescopes have then allowed determining the three dimensional, 
Keplerian orbits of stars within ten light hours of the compact radio source at 
the center of the Milky Way (see Figure [J). A few of them have already passed 
the pericenter allo wing iGenzel and Karad (120071 . and references therein) and 
Ghez et al.l (120031 ) to measure a central mass of ~ 3 x IO^Mq. 
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Figure 1: Left panel: a diffraction limited image of Sagittarius A* (~ 0.05 arcsec resolution) 
from the 8m ESO Very Large Telescope, taken with the NACO AO-camera and an infrared 
wavefront sensor at 1.6/2.2/3.7 /i m. The central black hole is located in the centre 
of the box. NACO is a collaboration between ONERA (Paris), Observatoire de Paris, 
Observatoire Grenoble, MPE (Garching), and MPIA (Heidelberg). Right panel: positions 
on the sky as a function of time for the central stars orbiting the compact radio source 
Sagittarius A*. T he data are from the UC LA group working with the Keck telescope. 
Figures taken from Genzel and Kara^ ( 2007 ) and references therein. 



BHs are currently measured in active galaxies too. The first attempts 
used direct spectral fitting of the Big Blue Bump, the optical/UV/soft X- 
ray continuum thermal emission from an optically thick, geometrically think 
accretion disk. More secure mass determinations are currently based on 
the Reverberation Mapp ing technique and extrapolations of the latter (see 
Ferrarese and Fordll2005l for a full description of such methods). 

Given the overwhelming evidence for a strong connection between BHs 
and their hosts, in particular all massive galaxies have a central massive 
BH, most theoretical models of galaxy fornaation now include recipes to 



grow BHs with i n galaxies (e.g.. iGranatq et al.l. l2004l: iDi Matteo et al.l. 12005 : 



Vittorini et al.. 2005: Bower et al., 2006: Croton et al.. 2006: Hopkins et al 



2006a 



Lapi et al 



2004, 2006, 2008 



20061: iMalbon et al.l . 120071 : Monaco et all . 120071 : iMenci et al. 



Cook et all . l2009f ). However, as explained below models 
differ significantly in the physical recipes adopted to evolve BHs, and al- 
though sometimes the observables reproduced are the same, the underlying 
physics might be widely different. 

It is therefore of extreme importance to probe BH evolution in ways 
independent of specific models. In this Review I will show how along the last 
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years, several groups have been able to develop and refine semi-empirical 
techniques to constrain BH evolution from first principles. The Review will 
focus on the latest statistical results on the BH/AGN population and show 
how these information can be interconnected to derive interesting and firm 
constraints on how BHs must have evolved along cosmic time. 

This Review will therefore r iot discuss details of BH/AGN emission pro- 
cesses (see the semi nal paper by ReesI 1984 ) , nor t he problem of seeding galax- 



ies with BHs (e.g., 
unification models 



Haiman and QuataertI 1200411. nor the problem of AGN 



e.g., lUrry and Padovanl Il995l : iMaraschil . 120021 : lElitzurl . 



20081 . and references therein). On the other hand, it will cover aspects of 



BH evolution which have not been properly discussed before in a coherent 
framework. 

In § [2] all the main correlations which have been measured between BH 
mass and host galaxy properties will be discussed, putting more emphasis 
on the most recent empirical results still not covered by previous reviews on 
the subject. In § [3] it will be shown how to measure the comoving number 
of BHs, and §§IHE1 and [H] will discuss why the latter information represents 
an invaluable tool to understand how BHs have grown to the population we 
observe today. §[7] will instead discuss how quasar clustering provides many, 
useful and independent constrains on BH evolution. A brief overview of how 
current models of galaxy formation are currently treating BH evolution, will 
be given in § [SI This aspect has not yet been systematically addressed by 
previous reviews on BHs, and it is useful to show how current open problems 
in galaxy-BH co-evolution models can be naturally and simply addressed via 
the techniques outlined in this Review. The main conclusions of this whole 
discussion will be listed in § [91 



2. The tight relation between Black Holes and their host galaxies 
in the local Universe 

The presence of nuclear BHs in the majority of galaxies (or at least 
bulges) has become an accepted paradigm. The nuclear kinematics of all 
nearby galaxies for which data of the necessary depth and spatial resolu- 
tion have been collected, show the clear signature of central mass concen- 
trations beyond what can be attributed to known stellar populations. The 
implied BH masses Mbh are found to correlate with several global properties 
of the host galaxies, including the mass, Mbuige (or luminosity Lbuige), veloc- 
ity dispersion, a, and light concentration Cbuige of the spheroidal component 
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surrounding dark matter halo (e.g., Ferraresj . 2002a ). In view of the clear 



,e.g., iFerrarese and Merrittl. l2000l: iGebhardt and et aU l2000l : I Graham et al.l . 
2OO1I iMarconi and Huntl . l2003l : iHaring and RiX l2004r) . and the mass of the 



importance of these correlations, it is worth discussing them below in some 
detail before progressing further in the Review. 

2.1. The Correlation between Black Hole mass and the luminosity /mass of 
the bulge 



Kormendy and Richstond (119951 ) reviewed the status of BH demography 



up to that time. Using a sample of only 8 well measured BH mass mea- 
surements they pointed out a tight c orrelation, with the BH mass being 



2-3% of the stellar mass of the host. iMagorrian et al.l (119981 ) constructed 
2D dynamical models for a sample of 36 nearby galaxies finding that the 
models required at the 95% confidence level the mass of the BH to be lin- 
early correlated with the mass of their hosts as Mhh^ O.OOGMstar- Updating 



the catalog of BH detections, Kormendy et al. (2000) concluded that the 



BH mass correlates much better with bulge luminosity than with the total 
luminosity of the galaxy, and that BH mass correlates with the luminos- 
ity of the high-density central component in disk galaxies independent of 
whether that component is a classical bulge (essentially a mini-elliptical) 
or a " pseudobulge" (believ e d to form via inward transport of disk mate- 
rial). iMcLure and DunlopI (120011 ) were the first to point out that if the 
sample is restricted to elliptical galaxies, the scatter is further significantly 
reduc ed. Despite the strong improvernents in bulge-bar-d i sk dec ompositions 



e.g., iPeng et al.l . l2002l : iGadottil . l2008l : iLaurikainen et all . 120091 ). due to the 



still small sample of local galaxies with secure BH mass measurements, it is 
still not clear whether the reduced scatter in ellipticals is mainly due to uncer- 
tainties associated in determining bulge masses in later-type galaxies, or t o 
an intrinsic physical effect (see also discussions in IFerrarese and Fordll2005l ). 
McLure and Dunlop (l200lh also showed that both active and inactive galax- 
ies follow similar correlations both in slope and normalization, although the 
int rinsic scatter around the correlations for inactive galaxies is much lower. 



Marconi and Huntl (120031 ) used K-hand images from the 2MASS database. 



and performed an accurate bulge/disk decomposition for the nearby quies- 
cent spiral galaxies with dynamically measured masses. They confirmed a 
decrease of the intrinsic scatter around the correlations when only early- 
type galaxies were considered. However, their total i^'-band magnitudes 
were derived by extrapolating the image profiles to very large distances on 
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implyinK an average in crease of brightness of about 0.5 magnitudes (e.g. 



Batcheldor et al. 



2004). 



Graham! (120071 ) critically revisited all the differing Mbh--^^ relations dis- 
cussed above by McLure & Dunlop, Marconi & Hunt and others. He con- 
sidered several previously not properly addressed issues finding that, when 
these adjustments are consistently applied to all the previous samples, results 
in relations which predict similar BH masses. His final result in the K-hand 
reads 

log(Mbh/MQ) = (-0.37 ± OM){Mk + 24) + 8.29 ± 0.08 , (1) 
with an intrinsic scatter of 0.33 dex. 



Gultekin et al.l (120091 ) recently recomputed the relations between BH and 



host-galaxy bulge velocity dispersion and luminosity, using the updated sam- 
ple of 49 Mbh measurements and also adopting 19 upper limits. Restricting 
to early-type galaxies, they find 



log(Mbh/Mo) = (8.95 ± 0.11) + (1.11 ± 0.18) x log(Lv/10"L£ 



Q,V) ) 



(2) 



with an intrinsic scatter of 0.38 ± 0.09 dex. The Authors claim their results 
to be insensitive to a wide range of assumptions about the measurement 
uncertainties and the distribution of intrinsic scatter. They also discuss how 
restricting the sample to the objects with resolution of the BH's sphere of 
influe nce might bias the results , a conclusion at variance w ith previous studies 



e. g.. iMarconi and Huntl. l2003l : iFerrarese and Fordl . l2005l ) 



Haring and Rixl (120041 ) approached the correlation between BH mass and 
host stellar mass in a way independent of photometric measurements, i.e., 
utilizing a total of 30 galaxies with bulge masses derived through Jeans 
modeling or adopted from djTiamical models in the literature. They found 
Mhh~ 1.4 X 10~^Mbuige with an observed scatter of < 0.30 dex, a fraction 



of which attributed to measurement errors. Their results confirm and refine 



the previous work of iMarconi and Huntl (120031 ). 

In sum, all of the most recent and detailed investigations have settled 
to define a closely linear correlation between BH mass and host galaxy 
mass/luminosity over several decades, with a factor of proportionality of 
2 X 10-3, e.g., Mbh^ 2 X lO'^Mh 



' bulge ■ 



2.2. The black hole mass velocity dispersion correlation 

The tight correlation between BH mass and the velocity dispersion a 
of the host galaxy is probably the correlation which received most atten- 



6 



tion in the literature as it provides a direct proof of the strong hnk be- 
tw een BH mass and potential well of the galaxy. I t was fi rst pointed out 
bvlFerrarese and Merritt fl2000h : iGebhardt and et al.l fl2000f ): iTremaine et al 



J2OO2I) , although with some discrepancies in the basic resu lts fsee alsolMarconi et al, 
2004) . and we refer to those papers and to the review by iFerrarese and Ford 
(I2OO5I ) for extensive discussions on this topic. 



One of the latest determinations of the Mbh-o" obtained by iTundo et al 



(120071 ) yields the result 



log(Mbh/M0) = (8.21 ± 0.06) + (3.83 ± 0.21) log (a/200 km s"^) , (3) 

with an intrinsic scatter of 0.22 ± 0.06 dex. These Authors also point out 
that the a-L relation of local early-type galaxies, including those having BH 
mass measurements, is different from the average a-L relation of early type 
galaxies calibrated in the Sloan Digital Sky Survey (SPSS). Thi s difference, 
they conclude, might bias the local BHMF. iLauer et al.l (l2007al ) also stress 
that the curvature in the a-L relation for massive early-type galaxies implies 
that the Mhh - a rela tion is unlikely to be a single power law. On similar 
lines, IWyithd (120061 ) suggested through a Bayesian analysis that a quadratic 
fit is a better description of th e co r relati on between BH mass and a in the 
Tremaine et al. (2002) sample. IHuI (120081 ) reinvestigated in detail the Mbh-cr 
relation, finding similar results to Eq. ([3]) for the ellipticals, but also pointing 
out that pseudo-bulges and core early-type galaxies seem to show significant 
departures from it, a possible signature of their different formation histories. 

2. 3. Redshift Evolution of the correlations 

Both the local characterization, and the cosmic evolution of the above 
mentioned scaling relations are subject to intense observat i onal studie s (e.g 



Shields et al.l. 20061: iPeng et al.l . l2006l : iLauer et al.l . l2007bl : iTreu et al.l . 12007 : 



Shankar et al.l . l2009al ). as they provide powerful cons traints on theoretical 
models for the co-eyolution of galaxies and BHs (e.g.. iGranato et al. 
Vittorini et aD.l2005l:lHoDkins et al.l. l2006al:ICrotonet al.l. l2006[ 



20061: ICavaliere and Mencil l2007l: iMonaco et all . l2007l : iMaruUi et all . l2008h . 



Treu et al.l (120071 ) and IWoo et al.l (120081 ) have randomly compiled from 



2004 



De Lucia et al 



the SDSS Data Release 4 a sample o f abou t 20 Seyfert galaxies in the redshift 
range 0.37 ^ z < 0.57. IShen et al.l (120081 ) estimated the Mbh-c relation for 
a larger sample of active galaxies up to 2; = 0.452. While the latter claim 
that no significant evolution in the Mbh-o" relation is detectable from their 
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sample, IWoo et al.l (120081 ) confirm the results by iTreu et al.l (120071 ) that a 
significant increase of ~ 0.2 dex in BH mass at fixed velocity dispersion must 
occur within z = and z ~ 0.5. However, systematic uncertaint i es ma y 
affect these estimates; for example, as also discussed by IWoo et al.l (120081 ). 



especially in galaxies with lower BH mass the host galaxy contribution to 
the 5100 A luminosit y may lead to an overestimation of the true BH mass. 



Peng et al.l (120061 ) collected a sample of 31 lensed and 18 unlensed AGNs 



at redshifts z > 1.7. They measured rest-frame R-band luminosities from 
H-band fluxes and BH masses by applying virial relations based on emission 
line widths. They found that the BH-to-host galaxy luminosity at 2; ~ 2 
is about the same as that at 2; ~ 0. Therefore, once the observed rest- 
frame luminosity has dimmed through substantial evolution to z ~ 0, at 
fixed BH mass the ratio BH-to-host luminosity grows significantly, and the 
resulting BH-luminosity normalizat ion is several t i mes h igher than the local 
one. Similar results were derived by lMcLure et al.l (120061 ) . who measured the 
BH-to-host galaxy mass ratio in a sample of radi o-loud AGNs in the r edshift 
range < z < 2 finding Mbh/Mstar oc {1 + zf. IShields et all (l2006h found 



that the CO emission lines in a sample of z > 3 quasars are very narrow, 
suggesting bulge masses about an order of mag nitude lower than rn easured 
in the local Universe, a t fixed BH mass (see also lCoppin et al.ll2008l ). 



Lauer et al.l (l2007bl ) have discussed several possible biases which may 
seriously affect these findings. At high redshifts a sample will be biased 
toward the most luminous AGNs and more massive BHs. Given the observed 
scatter in the local relations, especially significant in the BHhost luminosity 
relation, these massive BHs will be preferentially associated with the less 
massive but more numerous galaxies yielding a false signal of evolution. 

Several approaches have been pursued in the literature to address the is- 
sue of redshift evolution in the relations adopting methodolog ies which min- 
imize the impact of observational biases. iMerloni et al.l (120041 ) compared the 
accreted BH mass density in AGNs with the co smological global star forma- 
tion rate density (see also iHaiman et al.l 120041 ). Although their conclusions 
depend on additional assumptions about the fraction of the star forming 
galaxies which are linked to BH growth at a given redshift, they were able 
to rule out a strong evoluti on at a high c o nfidenc e level; we will further 



discuss their results below. iHopkins et al.l (l2006d ) also found, through a 



model-independent integral constraint based on comparing the AGN density 
with the luminosity and mass functions in different bands within the redshift 
range z = 2, very significant evidence for an extremely mild evolution. 
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Shankar et al.l (j2009al ) utilized the local velocity dispersion function (VDF) 
of spheroids, together with their inferred age distributions, to predict the 
VDF at higher redshifts {0 < z < 6). They then compared via a redshift- 
dependent Mbh-o" relation the cumulative BH mass density inferred from the 
AGN luminosity function with that associated with the VDF at each red- 
shift. This comparison, insensitive to the assumed duty cycle or Eddington 
ratio of quasar activity, yielded a very mild redshift evolution, in agreement 
with previous studies. 

2.4- Other correlations between Black Holes and host galaxies 

With a more numerous sample of galaxies with updated estini a tes o f 



galaxy distances, BH masses, and Sersic indices n, iGraham et al.l (120071 ) 



found evidence for a quadratic Mbh-n relation, of the type 

log(Mbh/MQ) = (7.98±0.09) + (3.70±0.46) log(n/3)-(3.10±0.84)[log(n/3)]^ , 

(4) 

with an intrinsic scatter of 0.18!!^ogg dex and a total absolute scatter of 0.31 
dex. When the quadratic relation is extrapolated, it predicts BHs with 
masses of ~ 10'^ Mq in n = 0.5 dwarf elliptical galaxies, and an upper bound 
on the largest BH masses in the local Universe equal to ~ 1.2 x 10^ Mq. This 
relation when adopted to compute the statistics of local BHs, yields quite 
dif ferent results with r e spect to other methods, as discussed in § [31 



Marconi and HuntI (120031 ) explored the addition to the Mbh-cr relation 



of the parameter (the effective radius) to create a "fundamental plane" 
for BHs. From their 27 "Group 1" galaxies, they constructed a relation 
between Mbh and and reO"^ (proportional to the virial bulge mass), which 
resulted in an intrinsic dispersion of ~ 0.25 dex, comparable to the one in 
the Mbh-(T rela tion. Allowing t he exp onents on the and a terms to vary 
inde pendentlv. iHopkins et al.l (l2007al ) used the same 27 Group 1 galaxies 
from iMarconi and HuntI (120031 ) along with some updated measurements, to 
find the best-fit relation 



log(Mbh/M( 



(8.33 ± 0.06) + (0.43 ± 0.19) log[re/3 kpc] 
+ (3.00 ± 0.30) log[a/200 km s"^] , 



(5) 



with an intrinsic scatter of 0.21 dex and a total scatter of 0.30 dex (see 



Graham 2008). 



Barway and Kembhavil (l2007l ) made the interesting claim that a combi- 



nation of two photometric parameters, namely the effective radius r^. and the 
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mean effective surface brightness /ig = — 2.51og/e, can be used to predict 
BH masses with a higher degree of accuracy than single quantities such as 
luminosity or velocity dispers ion. The optim al (B-band) solution using all 
18 data points as reported by iGrahamI (120081 ) is 



log(Mbh/Me 



&) 



(8.18 ± 0.09) + (3.15 ± 0.33) log[re/3 kpc] 
-(0.900.18)[/ie,s-21.0], 



(6) 



with a total scatter of 0.32 dex. 

The hypothesis of a BHFP which might be more fundamental than the 
Mbh-o" relation is obviously of great interest. However, later studies have 
questioned the actual existence of a BHFP. Graham (2008) showed that the 
barred early-type galaxies contained in the sample systematically deviate 
from the Mbh-cr relation, causing the residuals about the Mbh-o" relation to 
co rrelate with r^. giving ris e to the apparent BHFP (s ee, however, discussion 



Younger et al.l l2008l ). lAUer and Richstond (120071 ) investigate a variety 



m 

of correlations between bulge properties and BH mass, concluding that the 
strongest correlation exists between the bulge gravitational binding energy 
Eg, as traced by the stellar light profile, and the BH mass, of the type 
Mbhoc Eg'^. They, however, do not find evidence for a significant correlation 
between between BH mass, a, and Ve in early-type galaxies, with their best-fit 
rel ation being MbhCx: 



,3.73„ 0.05 
' p. 



Kormendy and Bended (120091 ) have adopted accurate measurements of 



the stellar "light deficit" Ldef associated with the cores of a sample of el- 
liptical galaxies characterized by a surface brightness covering a quite large 
dynamic rang e. They find, in agreement with several previous studies (e.g., 
Merrittll2006l . and references therein), that L^ei correlates well with the BH 
mass, with a similar, or even smaller scatter than the one in the Mbh-o" 
relation. They claim that this might be associated to dry mergers events, 
during which energy can be efficiently transferred to the central stars via the 
anisotropi c gravitationa l radia tion emitted after the coalesce o f the binary 



BHs ( e.g., ISesana et al.l . l2009l . and references therein). Finally, ISeigar et al. 



(I2OO8I ) claimed another correlation between the spiral arm pitch angle (a 
measure of the tightness of spiral structure) and the BH mass in the nuclei 
of disk galaxies. 



10 



3. The local mass function of Black Holes 



It is now widely accepted that knowing the exact number density of galax- 
ies in all bins of stellar mass (or other property) is fundamental to discern 
among successful models of galaxy formation. For example (see § E]), most 
semianalytic models (SAMs) of galaxy formation now agree that some sort of 
strong stellar and/or AGN feedback must have been responsible for shaping 
the exponential fall-off of the stellar mass function at high stellar masses. 

For similar reasons, it is essential to define the BH mass function (BHMF), 
i.e., the comoving number density of BHs per bin of mass in the local Uni- 
verse. Any model of galaxy formation and evolution which attempts to pro- 
vide a complete picture of galaxy evolution, must also address the problem 
of birth and growth of the massive BHs in galaxies, explain the nature of 
their tight relationship with their hosts, and predict their correct number 
and mass at z = 0. 



The standard procedure (first pioneered by lSalucci et al.lll999l and lFerrarese 



2002bj) in calculating the BHMF assumes that all galaxies with a signifi- 
cant bulge component host a BH, the mass of which can be inferred from 
the BH— galaxy scaling relations discussed in § |2j Knowledge of a local 
galaxy luminosity, velocity dispersion or light concentration function, com- 
bined with the corresponding BH scaling relation, will then lead to a BHMF 
through a sirnple change of variables. The first derivations of the BHMF (e.g., 



Salucci et al 



1999; Ferrarese. 2002b: AUer and Richstone. 2002: McLure and Dunlop 



2004 : lMarconi et all . 12004 : IShankar et all . 120041 : iTundo et all . 120071 : iGraham et al.l . 
20071 ) combined the galaxy velocity dispersion function (VDF) with the Mbh- 
a relation, and the spheroid luminosity function with the Mhh-MhuigR re lation. 



Here it is provided a brief compendium (see lShankar et al.ll2009bl ) of the 
BHMFs derived from applying different BH mass estimators. The aim here 
is to provide the reader with a compact overview of the most important steps 
required to generate a BHMF. 

The BHMF is computed by converting the galaxy distribution $(|/), 
which expresses the measured galaxy number density as a function of a given 
variable y (e.g., the stellar velocity dispersion or bulge luminosity), into a 
BHMF by adopting the corresponding empirically measured Mbh-?/ relation. 
The conversion is then performed via the convolution 



$(M, 



bhj 



^/2^TT]■ 



exp 



(Mbh - [g + by]y 
2ry2 



dy 



(7) 
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which takes into account the intri nsic scatter n in the M - h^-y relation, de rived 



from the direct data fitting (e.g., iTremaine et aLl . l2002l : iGrahaml . 120071 ). 



The BHMF can be derived by couphng the statistical information on 
local luminosity functions of galaxies with relationships among luminosity (or 
related quantit ies such as stellar m ass and velocity dispersion) and the central 



BH mass (see ISalucci et al.lll999l ). Because the BH mass correlates with 



the luminosity and velocity dispersion of the bulge stellar population, it is 
necessary to first separate the luminosity functions for different morphological 
types, that have different bulge-to-total luminosity ratios. It is convenient 
to use galaxy luminosity functions derived in bands as red as possible, where 
the old bulge stellar populations are more prominent. 

The spheroidal luminosity function $(Lbuige) for lenticular galaxies can 
be computed from the luminosity function of early- type galaxies, $e(-^), as 



$(L 



bulge) = T^^Mb X L) , (8) 
Je + /so 

f Yu and Tremaine 2002h. where fp; an d fso (which can be taken, for example, 
from Table 1 of iFukugita et al.lll998l ) are the numerical fractions of ellipti- 
cal and lenticular galaxies, respectively, to the total galaxy population, and 
b = Lbuige/-^, is the luminous bulge-to-total ratio for 5*0 galaxies. The contri- 
bution of ellipticals to the spheroidal LF can be computed by simply setting 
5 = 1 and replacing /so with /e in equation ([H]). Similarly, one can compute 
the luminosity function of bulges of spirals of different Hubble type by re- 
placing ^e{L) in equation ([8]) with the late-type galaxy luminosity function, 
and plugging in the appropriate /-weights for later type gala xies. 



A somewhat different approach was recently adopted by iGraham et al 
( 120071 ). First a BH mass via Eq. (jl]) and the appropriate statisti cal weight 
were associated to each spheroid in the Millennium Galaxy Catalog (jPriver et al 
20071 ). The latter was computed as the space density of of the appropriate 
spheroid type in the specified color and luminosity intervals, divided by the 
number of galaxies which contributed to that interval. The BHMF was then 
computed separately for early and late-type galaxies by summing the distri- 
bution of appropriate BH masses times their weights (see lGraham et al.ll2007 
for further details). 

When computing the BHMF from the velocity dispersion function of 
galaxies, instead, no correction for bulges is usually applied, at least for 
the bulge-dominated galaxies where the velocity dispersion a should be a 
good indicator of the true a of the bulge component (see discussion in 
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Ferrarese and Ford l2005l. and references therein) . Several grou ps therefore 



(e.g., iMarconi et all |2004J : IShankar et al.l . 12004 ). adopted the ISheth et al. 



(I2OO3I ) early-type velocit y dispersion function, directly converted into a BHMF 



via the Mbh-o" relation. IShankar et al.l (12004 ) also adopted the velocity dis- 



persion function computed via the bivariate equation 



(9) 



where pij is the fraction of sources in the lBernardi et al.l (120031 ) and lSheth et al 



(I2OO3I ) sample with velocity dispersion log cTj and luminosity Lj, normalized 



to the total number of sources with luminosity Lj, and ^{Lj) is the input 
luminosity function. The contribution of late-typ e galaxies t o the veloc- 



ity dispersion function is more difficult to assess. ISheth et al.l (120031 ) con- 
verted the late-type galaxy luminosity function into a velocity dispersion 
function by adopting a mean luminosity-velocity dispersion relation plus a 
m odel to correct f or gala xy orientation. The latter result was then adopted 
by lShankar et al.l (l2009bl ) to compute the contribution of the bulges of spirals 
to the BHMF. 

Figure [2] compares the local BHMFs obtained from using different rela- 
tions between BH mass and host galaxy properties. Each line in the Figure 
corresponds to a different BH niass est imator, as listed in the legend (we 
defer the reader to IShankar et al.ll2009bl for further details). 

The grey band of Figure [2] can be considered as representative of the 
mean and the systematic uncertainties of present estimates of the BHMF. 
The integrated mass density of the local BH population is pbh = (3.2 — 
5.5) X 10^ Mq Mpc^'^ (for h = 0.7). Figure [2] also presents two additional 
estimates of the BH MF. The dotted curve shows the estimate, for early- 
type galaxies only, by Hopkins et al.l ( 2007a ). based on the BHFP (Eq. [S]), 



which is in reasonable agreement with other determinations of the BHMF 
at the high mass end. The disagreement at the low mass end is due to the 
fact that the contribution from the bulges of spirals is missing from their 
estimate. Open circles show instead the estimate by I Graham et al. 1 (I2OO7I ) 
based on the correlation between BH mass and the Sersic index n (Eq. [1]). 
The latter includes a complete census of all galaxy types, and therefore the 
strong disagreement at the low mass end is a manifestation of a real mismatch 
in the results. This difference might be caused by both the use of a different 
relation (Eq. [3]), and, as explained above, a different method with respect 
to that outlined in Eq. ([H]), although this is difficult to posteriori. 
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Figure 2: Comparison among estimates of the local black hole mass function taken from 
Shankar et alj (|2009br ). Lines show estimates using various calibrations of the Mb 



Mbh - cr, or M, 



bh 



A^buige relations, assuming a 0.3-dex intrinsic scatter in all cases. 
The grey band encompasses the range of these estimates. The thick dotted line shows 
the determination of the early- type black hole mass function by iHopkins et al. ( 2007a ) 
using the blac k hole "fundamental plane" , and open circles show the determination of 



Graham et al.l (j2007[ ) using the relation between black hole mass and Sersic index. 
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Nevertheless, it is interesting to note that all estimates agree in providing a 
similar peak in the BH mass distributions, which is important to constrain 
BH accretion models, as we will discuss below. 

3.1. Black Holes in Later type galaxies 

One important caveat in determining the overall census of local BHs, is 
the actual presence of BHs in bulgeless galaxies. In the previous section, 
the assumption has been made that all galaxies with a bulge component 
possess a central BH. However, there is no firm observational evidence to 
support this claim. On one side, bulgeless AGNs are not uncommon (e.g., 
Gebhardt et al.l . 1200 ll : ISatyapal et al.l . l2008l ) , suggesting that current BHMF 
estimates might be missing the contribution of a BH population in late type 
spiral galaxies (although, based on current estimates, these BHs are 1 — 3 
orders of magnitude less massive than those discussed in this Review). On 
the other, the ubiquitous presence of BHs, especially in low luminosity dwarf 
ellipticals, has not been proven: for ins tance, it is questio nable whether a 
BH exists in the E5 galaxy NGC 205 fIValluri et al.ll2005f ). A number of 
recent studies, mostly with the Hubble Space Telescope, have shown that the 
majority of galaxi es, regardless of Hubble type, host morphologically dis tinct 



stellar nuclei (e.g.. lCote et al.l . l2006l : iFerrarese et al.l . l2006l : lBoekeii . 120071 ). the 
only exception (as a class) being bright galaxies with Mb ^ —20. 5, although 



Seth et al. 



cases of galaxies with both stellar nuclei and BHs do exist (e.g 
2008h . 

Evidence fo r BHs in snaaller g alaxies also come from the Chandra X-ray 
observations of I G alio et al.l (120081 ) . who showed that a large fraction of early- 
type galaxies in the Virgo Cluster (AMUSE- Virgo) have a significant X-ray 
emission, consistent with what expected from accretion onto a central BH. 
Satyapal et al.l (120081 ) reported the discovery of AGN activity in several local 
spiral galaxies, and concluded from their estimation of bolometric luminosity, 
that the average mass for the BHs residing in thes e galaxies should be around 
4 X IO^Mq, if they emit at the Eddington limit. IChosh et aD (j2008l ) found 
secure X-ray AGN detections in several nearby, face-on spiral galaxies. 

It remains still open the major challenge of understanding how BHs have 
grown and whether they followed the assembly history of their host galaxies 
(see § [8]). It is discussed in the following sections how, in ways indepen- 
dent of specific models, important constraints on how BHs have actually 
evolved within their host galaxies and dark matter halos can be derived by 
using global information on the statistics of local BHs, the statistics of active 
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galaxies and their correlated clustering strengths at different luminosities and 
redshifts. 



4. The integrated mass density 



The classic ISoltanI (119821 ) argument relates the integrated BH density to 
the integrated emissivity of the AGN population, setting interesting con- 
straints on the average radiative efficiency e of BHs and its possible de- 
pendence on mass and/or redshift. If the average efficiency of converting 
accreted mass into bolometric luminosity is e = L/MinflowC^, where Minfiow 
is the mass accretion rate, then the actual accretion onto the central BH is 
-^bh = (1 ~ e)^infiow, where the factor 1 — e accounts for the fraction of the 
incoming mass that is r adiated away instead of being added to the BH (e.g., 
Yu and Tremaind 120021 ). The rate at which mass is added to the BH mass 
function is then given by 



dp] 



bh 



1 



dt 



^{L)Ld\ogL, 



(10) 



with $(-L) the bolometric AGN luminosity function. 

Here I show a simple plot which conveys the basic idea of ISoltanI (Il982l )'s 
argument, ^ yhile I refer to o ther papers for a mo r e detailed discussion on the 
topic (e.g., iMarconi et all . I2004J : iMerloni et all . I2004J : IShankar et all . 12004 



20081: 



Tamura et al.l.l2006l: Hopkins et all . l2007bl : iMerloni and Heinzl . l2008l : [Silverman et al.l . 



Shankar et allJioogbh . Figure [3] shows the integrated accreted mass 
on BHs as a function of redshift following Eq. (|T0|) . for different values 
of the radiative efficiency e, as labeled. The grey stripe at z ~ in- 
dicates the systematic uncertainties in determining the local BHMF (see 
Figure [2]). It is clear from the Figure that the match between the local 
and accreted mass densities constrains the mean radiative efficiency to be 
e < 0.1. However, s till significant uncertainties on the bo l ometric correc- 



the local BHMF on the other (e.g., IShankar et all . l2009bl . and references 



tion on one side (e.g..lElvis et all. Il994l: iMarconi et al.l . |2004J : iHopkins et al. 
2007b : Vasudevan and Fabian . 20071 . 2009), and syst e matic uncertainties on 



therein), prevent a more precise comput atio n of the radiative efficiency (see 
also discussions in lShankar et al.ll2009bl an d lYu and Lull2008l ). 

Extended versions of the ISoltanI (Il982l ) argument have been successfully 
applied to constrain other aspects of BH evolution. An upp er limit to the 



degree of evolution of the Mbh-o" relation has been set by IShankar et al 
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Figure 3: Left panel: integrated accreted mass on black holes as a function of redshift 
following Eq. PU)) . for different values of the radiative efficiency e, as labeled. The grey 
stripe at z ~ indicates the systematic uncertainties in determining the local BH mass 
function (see Figure [2]). The match between the local and accreted mass densities provides 
interesting constraints on the mean radiative efficiency of BHs. Right panel: evolution in 
the best-fitting stellar mass density as a function of redshift derived from the mean BH 
accretion history rescaled by a r edshift-depend ent constant {solid line), c ompared with 
differe nt data sets as collected in iMerloni et al. Figure taken from iMerloni et al 

t004) . 



( l2009al ) by comparing the BH mass function associated with the velocity 
dispersion function at each redshift Zc (corrected for galaxies with age z < Zc), 
and the accumulated total BH mass d ensity extracted frorn Eg. f[TO|): these 



Author s rea ch similar conc l usions as Merloni et al.l (120041 ) . iHopkins et al. 



(l2006c[ ). and IShankar et al.l (l2009bl ) who used somewhat different methods. 
Some of these studies have also shown that the inferred average history of BH 
growth apparently parallels the cosmological star formation rate, at least in 
spheroids, at nearly all cosmological times, suggesting that the two processes 
may well be intimately linked and that, on average, BHs and galaxies may 
coevol ve. The right panel of Figure [3] shows the result of IMerloni et al. 
( 120041 ). These Authors compared all the available observational data on the 
redshift evolution of the total stellar mass (data points in the Figure), with 
the integrated mass evolution of BHs, estimated from the AGN luminosity 
function (solid line). Assuming that the ratio of the BH mass density to the 
stellar mass in spheroids evolves as (1 + 2;)", they found that the match to 
the data requires a mild, although significant, evolution in the BH-spheroid 
mass ratio with a ~ 0.5 ± 0.2; this implies slightly larger BH masses at 
fixed host stellar mass at 2; > 0, in agreement with the evolution inferred 
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by IShankar et al.l (j2009al ) for the Mbh-o" relation. Also shown in the same 
Figure are the relative contributions from early-type (dashed line) and late- 
type galaxies (dotted line) , with their corre s pondin g 1-a confidence intervals. 



Merloni and Heinzl (120081 ). IShankar et al.l (l2008a( ). and ICattaneo and Best 



((20091) have also used the integrated emissivity of radio sources to constrain 
the kinetic efficiency of BHs, in the hypothesis that radio power is correlated 
with the kinetic power of the jets. 



5. The evolved Black Hole mass function 



Important additional insight into the overall evolution of BHs can be 
gained by comparing not only the total m ass densities, but also the shapes 



of the mass functions (jSalucci et al.lll999l ). In fact, the comparison between 



the predicted mass distribution implied by the AGN luminosity function 
(Eq. [H]), and the local BHMF can yield further interesting constraints on 
the mean Eddington ratio A at which BHs are accreting at. 

The predicted mass function n(Mbh,t) (in units of Mpc'^Mg^) implied 



nuitv eauation argument (e.s.. Cavaliere et al.. 


1971 


Small and Blandfordl. 


1992; 


Yu and Tremaine, 


2002; 


Marconi et al.. 


2004J: ^ 


tu and Lul. 


2004, 


2ooa 


Shankar et al.. 


2004; 


Hopkins et al.. 


2007b; 


Shankar et al. 


. 20091: 


) 



dn ,^ ^ 
^(M.„t) 



9((Mbh)n(Mbh,t)) 



bh 



where 



(Mbh) 



(11) 



(12) 



is the mean accretion rate (averaged over the active and inactive popula- 
tions) of BHs of mass Mbh at time t, and A = L/L^dd the Eddington ratio, 
with LEdd = 1-26 X lO^^(Mbh/M0) ergs"^ the Eddington luminosity (i.e., 
A oc L/Mbh). Eq. (fTTj) states that the average growth rate of all BHs is pro- 
portional to the function S'(Mbh, z, A), i.e., the fraction of BHs of mass Mbh 
active at redshift z and accreting at the Eddington rate A. This evolution 
is equivalent to the case in which every BH constantly grows at the mean 
accretion rate (Mbh). In practice, individual BHs turn on and off, and there 
may be a dispersion in Mbh values, but the mass function evolution depends 
only on the mean accretion rate as a function of mass. 
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Eq. (fTTj) can be further generalized to any distribution 



^(Mbh, z, A) = p(A, ^)A(Mbh, z) ; (13) 

this assumes all active BHs at redshift z to share the same mean Eddington 
ratio distribution p(A, z) (normalized to unity), and A(Mbh,2;) is the duty 
cycle, i.e., the total fraction of active BHs at redshift z and mass Mbh, ne- 
glecting for simplicity any further mass-dependence of the Eddington ratio 
distribution p(A, 2;). A physically consistent model must have A(Mbh, z^ <\ 
for all Mbh at all times. 

In models where all BHs grow with a single value of the Eddington ratio 
A, p(A, z) is just a Dirac delta-function p(A, z) = 5(A — A), while the duty 
cycle is simply the ratio of the luminosity and mass functions, 

where $bh(Mbh,2:) = n(Mbh, 2;)Mbh In(lO), and ^{L,z) is in units of Mpc"^ 
dex~^ 

Once an initial condition is specified, the solution to Eq. flTTl) can be 
performed either iteratively, updating the duty cycle computed from the 
AGN luminosity function at the appropriate timestep (Eq. [T3j) or in simple 



integral forms, at least in the case of single-A models (see, e.g., lYu and Lu 



2004 ). On adopting a broad input p(A, z') in which different fractions of BHs 
of the same mass are allowed to accrete at different values of the Eddington 
ratio A, Eq. fllip can be efficiently solved iteratively. At any given time, 
given the observed AGN luminosity function $(L, z\ the underlying BH mass 
function $bh(Mbh, 2:), and Eddington ratio distribution p(A, z^, the duty cycle 



A (M bh, z) is found fully specified by the equality (e.g.. lSteed and Weinbergl . 



20031 ) 



$(L,z) = y p(A,z) A(Mbh,2:)$(Mbh,2)rflogMbh, (15) 

given the input AGN luminosity function $(L,2;) and the Eddington ratio 
distribution p(A, 2). 

Figure H] shows the results from the first attempts to match the local 
BHMF with the BH mass function as predicted from Eq. ffTTj) with the as- 
sumption that all BHs accrete with a single val ue of the Edd i ngton ratio. 



The left panel shows the result worked out by ISalucci et al.l (119991 ). with 



the dotted and dashed lines referring to the predicted mass function for two 
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Figure 4: Left panel: predicted mass function from Eq. pT|) as derived from lSalucci et al 
(Il999l) assuming A = 1 [dotted line) and A = 0.2 (dashed line). The solid line is their 
estimate of the local black hole mass function derived from the Mbh-Mstar relation, while 
the filled symbols refer to their e s timat e of the local mass function derived from radio data. 
Figure taken from lSalucci et aL I (ll999l) . Right panel: cumulative blac k hole mass function 
{thick solid line) computed from Eq. (fTTj) bv lYu and Tremaind ( 2002 ) assuming A = 1 and 
e/(l — e) = 0.1 compared to local estimates of the local b lack hole mass function shown 
with dotted and dashed lines (see Yu and Tremaine 20021 for details). The thin solid line 
is their estimate of the predicted black hole mass function cutting out quasars radiating 
above a luminosity L = LEdd(10^ Mq). The lower panel shows the same quantities as in 
the upper panel normalized to unity. 



values of A = 1, 0.2, respectively, compared with their estimates of the local 
BHMF (solid l ine and filled symbols). T he right panel shows instead a similar 
comparison by lYu and Tremaind (120021 ). where the solid lines show the cumu- 
lative BH mass functions obtained assuming A = 1 and e/(l — e) = 0.1, com- 
pared to their estimates of the local BHMF shown with dott ed and dashed 
lines. These preliminary works found, in agreement with the ISoltanI (119821 ) 
argument discussed in § HI that the average BH radiative efficiency should 
be around e < 0.1, in the assumption that most of the BH mass in the local 
Universe was accreted during AGN visible phases. In addition, they were the 
first to point out that, once e is fixed using Eq. ([TU]) . the overall shape of the 
local BHMF can be naturally reproduced by direct integration of the AGN 
luminosity over time using Eq. (fTTj) , once the appropriate median A has been 
adopted. 
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The upper panels of Figure show the results of a reference model which 
evolves the predict ed mass function assum ing an Eddington ratio of A ~ 0.40 
and e = 0.065 (see Shankar et al. 2009bl for further details). The upper left 
panel of Figure [5] plots the duty cycle A(Mbh, z) as a function of mass for 
different redshifts, as labeled. The duty cycle for Mbh ~ 10^ is ~ 0.2 at 
z ~ 4 — 5, falls to 0.03-0.08 at 2; = 2 — 3 when quasar activity is at its peak, 
then drops to 0.003 at 2; = 1 and ~ 10"^ at z = 0. Below z ~ 3, the duty 
cycle rises towards low BH masses. This "downsizing" evolution, in which 
high mass BHs complete their growth early but low mass BHs continue to 
grow at late time s , is re quired by the observed luminosity function evolution 



e.g.. lUeda et al.l . l2003l ) in any model with approximately constant A. 



This model yields overall good agreement with the average determinations 
of the local BHMF. The upper right panel of Figure [5] plots Mbh$(Mbh), pro- 
portional to the fraction of BH mass per logarithmic interval of Mbh, which 
allows better visual comparison to the observed local BHMF and highlights 
the contribution of each BH mass bin to the total mass density at each time. 
The evolution of the mass function derived from the continuity equation is 
shown here at different redshifts from z = 3 (dot-dashed line) down to ^ = 
(solid line). Because of the luminosity-dependent density evolution in the ob- 
served luminosity function, the massive end of the BH mass function builds 
up early, and the lower mass regime grows at later redshifts. 

For Mbh > 10^'^ M0, the mass function is almost fully in place by 2; = 
1. At very high masses the model exceeds the observational estimate, but 
in this regime the estimate relies on extrapolation of the scaling relations, 
and it is sensitive to the assumed intrinsic scatter. In addition, the high- 
mass end of the predicted mass function is sensitive to the bright end of 
the AGN luminos ity function which might also suffer from so me systematic 



uncertainties (see IFontanot et al.l 120071 . IShankar et al.l l2009bl and references 
therein). The op en circles with error bars show the estimate of the local 
mass function by iGraham et al.l (120071 ). which cannot be reproduced even 
approximately with constant A-models. If this estimate is correct, then the 
low end of the luminosity function must be produced mainly by high mass 
black holes accreting at low A, so that the predicted growth of low mass BHs 
is reduced. 

Assuming that all BHs at all times radiate at a constant Eddington ratio is 
obviously an oversimplification. Some empirical evid ence for some sharp Ed- 
dington ratio distribution p(A, z) has been found by iKoUmeier et al.l (120061 ) 
who claim that most AGNs of different BH mass in the AGES survey at 
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Figure 5: Lower right panel: predicted duty cycle A(Mbh,2) on tlie basis of Eq. pT|) 
at different redshifts and BH mass, as labeled. Upper right panel: predicted black hole 
mass function ^{Mhh,z) from a model with a delta-function p{X). The result is shown 
at different redshifts, from z = 3 to 2: ~ 0, as labeled, in the plane $(Mbh, z)Mhh, which 
plots the mass density per bin of BH mass. The local mass function is shown with a grey 
band and has a shallower slope than the one predicted by this model at z ~ {solid line). 
The open circles with error bars show the estimate of the local mass function by Graham 
et al. (2007). Lower left panel: assumed Gaussian p(A) distribution with the peak steadily 
decreasing with decreasing redshift. Lower right panel: predicted black hole mass function 
<I>(Mbh, z) from a model with a p{X) as given in the lower left panel. Even in this case the 
result is shown at different redshifts, as labeled. 



0-5 ^ -2 ^ 4 all share a similar Gaussian-shaped p(A), constantly peaked 
around A ~ 0.25 and with a dispersion of ~ 0.3 dex. On the other hand, at 
least at 2; < 1, AGN Eddington ratio distributions have a wide distribution, 



with a large fraction of Seyferts r adiating at A < 0.1 (e.g. JHickox et al.l . 12009 ; 
Kauffmann and Heckmanl . l2008l ). Some studies also show evidence for strong 
redshift evolution and nailder mass- dependence of the typical Eddington ra- 



tio fe.g-.lHeckman et al 



Constantin and Vogeley 



2004: 



2006 



yestergaardl2004j:lMcLure and Dunlod. 12004 



Dasyra et al.l . 120071 : iNetzer and Trakhtenbrot 
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Figure 6: Left panel: average growth h istory of BHs with given starting mass at z = 3 
computed using the lUeda et al. luminosity function and e = 0.1, A = 1. The 

symbols, filled squares, empty squares, and stars, indicate when the black hole reaches 
90%, 50%, and 5% of its final mass, respectively. Right panel: average mean lifetime of 
active BH s (AGNs) as a fun ction of the relic BH mass at z = computed with e = 0.1, 
A = 1, the Ueda et al. ( 2003i) {solid line) and Bovle et al.l ( 2000l ) {dashed line) luminosity 
function. The dotted line is the corresponding Salpeter time. The scale of the y axis on 
the right is the average lifetime in units of the Salpeter time. The lines with the empty 
squares and triangles are the average lifetimes computed from the Ueda et al. (2003) 
luminosity function with e = 0.0 4, A = 0.1, and e = 0.15, A — 1.6, respectively. Figures 
taken from arconi et al.l (j2004l ). 



20071 : iRovilos and Georgantopouloa . 120071 : IShen et al.l . l2008l ) , though system- 



atic uncertainties in the reverber ation mapping techniques and extrapo lations 
of empirical virial relations (e.g.. iKaspi et al.l . l2000l : iBentz et al.l . |2006| ) make 
it hard to draw firm conclusions. 

The lower left panel of Figure [5] shows an input p{X, z) with negligible 
mass- dependence but a strong redshift evolution (as labeled in the panel) 
re quired to be simultaneo usly consistent with the high-A observed at 2; > 1 
bvlKoUmeier et al.l (120061) . and the low A ~ 0.01 observed in SDSS by, e.g., 
Kauffmann and HeckmanI (120081 ). At each redshift a significant spread in 



p(A, z) consistent with the observed one is also implemented. The predicted 
BH mass function $(Mbh, z) implied by this model is shown in the lower right 
panel of Figure O at different redshift from 2; = 3 to 2; ~ as labeled, plotted 
in the $(Mbh, z)M\^\^ plane. From the comparison with the local BHMF, it is 
clear that the latter model still predicts a similar break for the z = BH mass 
function as in constant A-models, because the average A is still about 0.3—0.5. 
However, much less numerous low mass BHs are produced when the median 
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A strongly decreases with redshift, because more of the low- 2; emissivity is 
mapped into more massive BHs {L oc AMbh), preventing the growth of the 
less massive ones. It is therefore clear that an extremely accurate estimate 
of the local BHMF is fundamental to discern among different BH accretion 
models (see § 18.21) . It is also possible to show that broadening the input 



p{X) distribution does not affect the downsizing of BHs at low redshifts, bu t 
produces an upsizing at higher redshifts (see also iMerloni and Heina l2008l ) . 
with the fraction of more massive active BHs significantly higher than low 
mass ones. 

Other interesting quantities which can be directly derived from the con- 
tinuity equation formalism developed above are the mean growth of BHs of 
a given mass at 2; = and the average visible lifetime of BHs. The former 
quantity can be directly drawn from direct integration on the mean accretion 
rate 

M^^{z)= / {M^^{z'))—dz' . (16) 



The re sult is shown in the left panel of Figure [6] as derived from lMarconi et al 
(I2OO4J ). w hich shows t h e aver age growth history of BHs from z = 3 computed 
using the lUeda et al.l (120031 ) luminosity function and e = 0.1, A = 1. The 
symbols, filled squares, empty squares, and stars, indicate when the BH 
reaches 90%, 50%, and 5% of its final mass, respectively. Here it is apparent 
the notion of downsizing, where the most massive BHs acquire most of their 
mass at 2; > 1, at variance with lower mass BHs. 

The corresponding time for which a BH of a given mass at 2; = was 
active in the past can be derived by direct integration of the duty cycle over 
time, as 



tvis[Mbh(2: = 0)] 



A(Mbh,t)cit. 



(17) 



The result, again taken from iMarconi et al. (l2004t ). is shown in the right 
panel of Figure [6] for a model with e = 0.1 and A = 1, u sing both the X- 
ray and opt i cal lu minosity functions of lUeda et al.l (120031 ) (solid line) and 
Boyle et al. ( 2000l ) (dashed line), respectively. Also the outputs of other 



models with slightly different accretion parameters are plotted in the same 
Figure for comparison. It is clear that lower mass BHs have longer visibility 
timescales than more massive BHs. Note, however, that the latter timescale 
strongly depend on the epoch of "first appearance" in which the BH of that 
final mass starts shining in the luminosity function, and therefore cannot 
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be precise estim a tes of BH lifetime s , as discussed by several Authors (e.g. , 
Hosokawa , 20021: Yu and Treraainj. 2002; Shankar et al. . 2004; Yu and Lu . 
2004J : iHopkins and Hernquistl . l2008h . 



6. Black Hole mergers 

So far we have assumed that BHs mainly grow through accretion. How- 
ever, observations show, and hierarchical galaxy formation models predict, 
that a significant fraction of galaxies experience mergers with comparably 
massive galaxies during their lifetime. At least some of these galaxy mergers 
are likely to be accompanied by mergers of the central BHs that they con- 
tain. In fact, as the galaxies merge, the BHs should sink toward the center 
of the new galaxy via dynamical friction where they form a binary, which 
can continue its decay by transferring angular momentum to the stars in- 
tersecting its orbit. If the binarys separation decreases to the point where 
the emission of gravitational waves becomes efficient at carrying away the 
last remaining angular momentum, the BHs coalesce rapidly. However, the 
mechanisms that shrink the orbits of the BHs with their small cross-sections 
to s uch small scales of the order of a ~ 1 parsec, are not fully understood 
(see lMerritt and Milosavljevial2005l for a review). 

In principle, including mergers in the overall evolution of the BH pop- 
ulation, requires inserting a mass and redshift-dependent source term in 
Eq. fill I) . During the evolution, mergers will then redistribute mass within 
the BH population, thus varying at any given epoch the predicted mass dis- 
tribution with respect to a pure accretion model. Moreover, mergers should 
not change t he integrated mass density, so they do not affect the integrated 



SoltanI (119821 ) argument discussed in the previous sections; in principle, grav- 



itational radiation during mergers can red uce the integrated mas s density. 



but it should not be a dominant effect (see lYu and Tremaindl2002i ). 

However, at present it is still difficult predicting precisely the BH mass 
function implied by gas accretion and mergers. On the one side, the actual 
rate of BH merging is still highly uncertain; on the other side, the indirect 
important effect that mergers might have on the mass- dependence of the 
radiative efficiency is still unknown. In fact, there is a tight connection 
between radiative efficiency and angular momentum of the BH. Most of the 
energy is radiated close to the innermost stable orbit, which ge ts closer to the 



centr al rotation axis if the BH is rapidly spinning (see, e.g., iBardeen et al 



19721 for the formalism). Therefore a high spin is paralleled by a high radiative 
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efficiency. BHs are believed to spin up or spin down after mergers and/or gas 
accretion, altliougli tlie details of such processes are too poorly understood 
to allow detailed modeling. The canonical value of the radiative efficiency 
competing to a non-rotating BH is e ~ 0.06, incr easing up to e ~ 30% for a 



rapidly spinning BH in equilibrium ( Thornel . ll974f) . The value of the radiative 
efficiency e < 0.1 inferred from the ISoltanI (119821 ) argument discussed in § H] 
(see Figure [3]), appear to indicate that on average BHs are not extreme 
rot ators, although the cav e ats d iscussed in § H] might bias this conclusion. 



Hughes and Blandfordl (120031 ) studied the remnants of merging BHs of dif- 



ferent sizes, spins and orbital parameters. They concluded that BHs are gen- 
erally spun down by mergers and only for a very narrow rang e of orbital pa- 



ramet ers merging is effective in rapidly spinning up the holes. iGammie et al. 



20041) carried out numerical simulations confirming the lHughes and Blandford 
(l2003l ) result that a BH is usually spun down by frequent minor mergers. 



wh ile major merge r s usu ally spin up the BHs. 



Volonteri et al.l (120051 ) studied the expected distribution of BH spins and 



its evolution with cosmic time in the context of hierarchical galaxy formation 
theories. They followed the merger history of dark matter halos as extracted 
from Monte Carlo realizations based on the Extended Press-Schechter for- 
malism, and followed the growth of BHs via gas accretion and mergers. They 
concluded that mergers do not lead to a systematic spin up or spin down of 
BHs with time, but actually the spin distribution retains memory of its ini- 
tial conditions. In their models, accretion is the main driver for efficiently 
spinning BHs up, and they find that if accretion takes place from a thin 
disk, then up to 70% of all B Hs can be maxi r nally rotating and have ra- 
diative efficiencies > 12%. The IVolonteri et al.l (120051 ) conclusion was based 
on the assumption that the BH spin aligns very rapidly with the accretion 
flow, although they also note that a randomly-or iented sequence o f accre- 
tion events is less efficient in spinning the BH up. iKing and Pringld (120061 ). 
strengthen the latter point, discussing through analytical methods how a ran- 
domly oriented sequence of accretion episodes can effectively keep the spin 
low, rather than alignin g the h ole very quickly and causing systematic spin 



up. iBerti and Volonteri] (120081 ) used results from numerical relativity simu- 



lations of BH mergers, coupled to the cosmological evolution of BH spin, to 
infer the evolution of the spin distribution. They conclude that mergers are 
very unlikely to yield substantial spins, unless alignment of the spins of the 
me rging holes wit h the or bital angular momentum is very effective. 



Shankar et al.l (l2009bl ) have illustrated the potential impact of mergers 
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on the predicted BH mass function using a simple, yet useful mathematical 
model that assumes constant probability of equal mass mergers per Hubble 



time, similar to the models of ISteed and WeinbergI (120031 ). Following the lat 



est results di scussed above on the re latively mild effect of mergers on the spin 



distribution, IShankar et al.l (l2009bl ) neglected the variations of radiative effi- 



ciency caused by the impact of mergers on the BH spin distribution, and also 



body interactions (e.g.. 


Huehes and Blandford 




2003; 


Gammie et al. 




2004; 


Islam et al.. 


2004; 


Yoo and Miralda-Escudel. 


2004 


Volonteri et al.. 2005 


Merritt and Milosavlievic 


2005; 


Yoo et al., ^ 


^007; J 


5esana et al.. 


2009). 


H of mass Mbh has a probability 




Shankar et al. 


(2009br) assumed that a B 



-Pmerg of merging with an equal mass black hole in the Hubble time tniz) 
(age of the Universe at redshift z). Therefore the fraction F of black holes 
that merge in a timestep At is given by 



F = P X 

merg ^ 



At 



(18) 



At each time ti, these Authors first advance the mass function to time t2 = 
ti + At with accretion only, then add to each bin of the mass function an 
increment (which may be positive or negative) 



A$(Mbh,t2 



Fx $(^,t2) 



-Fx $(Mbh,t2 



(19) 



where the second term represents BHs lost from the bin by merging and 
$ (^^^,^2) is calculated by interpolation. 

Figure [7] shows the evolution of the BH mass function resulting from a 

0.065 and A ~ 0.4. The net effect of merging is to 



model with P. 



merg 



0.5, e 



slightly lower the abundance of BHs below the peak of the mass function and 
to si gnificantly increase the abundance o f very massive BHs, as expected (see 
also lMalbon et al.ll2007l . lYoo et al.ll2007l ). Since a non-merger model already 
produces an excess of massive BHs relative to the local BHMF (see Figure [5]), 
adding mergers only makes the match to observations worse. However, the 
impact of mergers is evident mainly for Mbh > 10^ Mq, where the BHMF 
estimates are most sensitive to the adopted scatter in the black hole-host 
scaling relations and to the extrapolations of these relations up to the most 
luminous galaxies. Note that the value F mprg = 0.5 is an upp er limit for the 



merger rate of massive galaxies (see, e.g., iBundy et al.ll2009l and references 
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Figure 7: Evolution of the black hole mass function in a model with black hole mergers. 
Accretion-driven growth is computed assuming e = 0.065, A = 0.40, and each black 
hole has a probability Pmerg = 0.5 of merging with another black hole of equal mass 
per Hubble time tuiz). Squares show the z = Q predictions of the reference model (same 
accretion parameters, no mergi ng), and the long-dashe d line shows the z = mass function 
for a merger model with the Hopkins et al. ( 2007bh luminosity function and accretion 
parameters e = 0.09, A '--^ 1. 



28 



therein), and it was chosen to provide an upper hmit on the actual impact 
of mergers on the BH population. 

The impact of mergers on the predicted BH mass function is likely to 
be small compared to remaining uncertainties in accretion-driven growth, 
except perhaps for the rare, high mass BHs. The most interesting impact 
of BH mergers may arise indirectly, through their e ffect on BH spins an d 
thus on radiative efficiencies, as discussed above (e.g.. IVolonteri et all l2005l ). 
Assuming, for example, the radiative efficiency to increase with increasing 
BH mass, might be a viable solution to accommodate a sequence of numerous 
mergers with the extremely stee p fall-off of the local BHMF at the high-mass 
end (see also ICao and Lil 120081 ) . Potential evidences for some mass and/or 
redshift dependence of the radiative efficiency may already be available (see 

§[ZD- 

Figure [8] illustrates my summary of the main constraints which can be 
derived by comparing the predictions of accretion plus mergers models for 
BH evolution with the local BHMF. The overall B H mass densi ty con- 
strains the radiative efficiency as first pointed out by ISoltaru (119821 ). while 
the p eak of the mass dis tribution is able to constrain the mean Eddington 
ratio (jSalucci et al.l 119991 ). Matching the low end of the mass distributions 
can instead shed light on the alternative between models with low-mass BHs 
growing at a significant fraction of the Eddington limit, and models marked 
by higher- mass BHs growing at sub-Eddington regimes (see § Finally, 
the high-mass end, coupled with some information on the mass-dependence 
of the radiative efficiency, should provide useful constraints on the impact of 
mergers on the overall BH evolution. 



7. AGN CLUSTERING PROVIDES INDEPENDENT CONSTRAINTS 
ON BLACK HOLE EVOLUTION 

In the previous sections it has been shown that continuity equation ar- 
guments based on evolving the BH population linking BH growth to the 
observed AGN bolometric luminosity function yields interesting constraints 
when compared to the local BHMF. However, this technique can only con- 
strain mean quantities. For example, it cannot distinguish between a model 
with a constant A = 0.4 at all redshifts, or a model with Super-Eddington 
accretion (A > 1) at higher redshifts, steadily decreasing to sub-Eddington 
regimes at lower-z (lower panels of Figure [5]). On the other hand, there are 
no assumption-free direct observations of the radiative efficiency at high red- 
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Figure 8: Summary of the main constraints which can be derived by comparing the pre- 
dictions of accretion plus mergers models for black hole evolution with the local black 
hole mass function (grey area). The overall height under the curve constrains the ra- 
diative efficiency, while the peak of the mass distribution is able to constrain the mean 
Eddington ratio. Matching the low end of the mass distributions can instead shed light 
on the alternative between models with low-mass BHs growing at a significant fraction of 
the Eddington limit, and models marked by higher-mass BHs growing at sub-Eddington 
regimes. The high-mass end, coupled with some information on the mass-dependence of 
the radiative efficiency, should instead provide useful constraints on the impact of mergers 
on the overall BH evolution. 
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shifts (e.g.. lRafiee and Halll . l2009l : IWang et al.l . l2009l . and references therein), 
and a cons tant as well a s a redshift-dependent radiative efficiency can still 
satisfy the ISoltan argument. Breaking such degeneracies would re- 

quire knowledge of the BH mass function at higher redshifts, which is still 
not feasible given that it is still not clear whether the local relations between 
BHs and galaxies still hold at higher redshifts (e.g., IShankar et al.l . l2009al . 
and references therein). 

Quasar clustering provides additional, in dependent constraints on the 



main parameters regulating BH evolution (e.g.. lMo and Whitd . ll996l : lHaehnelt et al. 
19981 ) and their possible evolution with redshift and/or mass, thus help- 
ing in discriminating among the successful models. As first o utlined by 



Martini and Weinberg! ( 2001 ) and Hairnan &: Hui (2001; see also iLidz et al.l 



2006 1: [Hopkins et al. P2008': 'S hankar eraL 



20081 : iMarulh et al.l 12009. : Thacker et al 



2008b; White et al.ll2008: Bonohetal 



20091 ). the clustering is in fact a di- 



rect measure of the masses, and therefore number densities, of the halos 
hosting the quasars. In turn, the ratio between the quasar luminosity func- 
tion and the halo mass function provides information on the duty cycle, i.e., 
the fraction of halos which shine as quasars at a given luminosity, BH mass, 
or Eddington ratio. 

Many groups have now been able to measure the c lustering of AGNs at 



different luminosities, bands, scales, and redshifts (e.g.. lOverzier et al 



Grazian et al.l.l2004l:ICroom et al.l.l2005l:IConstantin arid Vogelevl. 



2003; 



2006; 



20041: Porciani and Norberel. 2006r da Angela et al 



Hennawi et al. 



20061: 



Porciani et al. 



Mvers et all l2007allb l: IPa dmanabhan et al.l l2008l: IShen et all. l2007l. 



2008 

Coil et all . l2007l . [2009; M ivaii et al.l . l2007l : iHickox et all . l2009l : IRoss et all . 



20091 . and references therein). AGN clustering is usually expressed in terms 
of the bias, i.e., the square root of the ratio between the matter and the mea- 
sured two-point correlation functions. It is beyond the scope of this Review 
to discuss and compare all the above mentioned empirical results. In the 
following we will just focus on the semi-empirical methods which have been 
adopted in the literature to actually use such measurements for deriving firm 
constraints on BH evolution. 

The classical modeling of quasar clustering by Haiman & Hui (2001) 
and Martini & Weinberg (2001) assumes a mean value for the duty cycle 
and derives the relation between quasar luminosity and host halo mass by 
monotonically matching their cumulative distribution functions. Formally, 
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this concept leads to 



^(x, z) dx 



-erfc 



In 



2^ Q?/min {•^min) 



10^^ J V21n(10)S, 



(20) 



here x = logL {L being the quasar luminosity), and y = log Mh {Mh being 
the halo mass). The quantity ^hiu Tz) is the comovin g number density of 



halos, in units of Mpc ^ dex (e.g., ISheth et al.l . 1200 ll ). while $(x, z) is the 



luminosity function of quasars in the same units of Mpc~^ dex~^. The quan- 
tity Ah{y,z) in Eq. fl20l) is the duty cycle, i.e., the fraction of halos which 
host quasars shining above a minimum luminosity x^i^ = log L^^^ at redshift 
z. Eq. fl20|) also takes into account a scatter of S (in dex) between quasar 
luminosity and halo mass. 

So the clustering clearly provides an independent way to constrain duty 
cycle Ah{y, z) of active halos and BHs, which can be directly compared with 
the A(Mbh,-2) predicted from continuity equation arguments (see Eq. [H]). 
Either method can be adopted independently to work simultaneous con- 
straints on the radiative efficiency, Eddington ratio distributions, and the 
underlying BH mass distribution at a given redshift; this ultimately provides 
in ways independent of specific models, powerful tools to dehneate a com- 
prehensive picture of BH evolution. Also, these determinations of the duty 
cycles can be further constrained by comparing with direct estimates of the 
fracti on of active galaxies iii complete samples above a certain luminosity 



e.g., iKauffmann et al.l . l2003l ). 



At each redshift Eq. (120|) defines the minimum halo mass ?/min correspond- 
ing to the minimum luminosity in the sample Xmin- The mean bias associated 
to the same subsample at redshift z is then 



{b) 



Cn dy^h{y,z)b{y,z) 



r, dy<l>h{y,z) 

^ ymin 



f211 



Martini and Weinberg! fl200lh and lHaiman and Hull fl200lh applied Eq. ([20]) 



showing that quasar clustering measurements can substantially narrow down 
the range of possible duty cycles. In fact, if quasars are long-lived they are 
rare phenomena that are highly biased with respect to the underlying dark 
matter; whilst if short lived, they reside in more typical halos that are less 
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Figure 9: Left panel: correlat ion length pred i cted fo r two different quasar lifetimes, as 
labeled. Figure taken from ,Haiman and Huil (j2001r ). Right panel: clustering length vs 
quasar lifetime tg and for different values of th e quasar number density <&, as labeled. 
Figure taken from 'Martini and Weinberg ( 200l[ ). The main point in both panels is that 
an higher quasar lifetime (i.e., an higher duty cycle) implies the quasars to be hosted in 
less numerous, more massive and more biased halos. 



strongly clustered. Figure [9] shows examples of their calculations in which 
the duty cycles are actually expressed in terms of the quasar optical lifetime 
tg, which is directly proportional to the duty cycle A^, for given approxima- 
tions as to the lifetime of the host halos. Both plots show that the clustering 
strength increases for increasing quasar lifetime tg and, at fixed tg, increases 
for decreasing quasar number density as labeled. This behavior is easily 
understood given that longer quasar lifetimes (i.e., an higher duty cycle), 
and/or lower quasar abundances, imply quasars to be hosted in less numer- 
ous, more massive, and more biased halos. 

Eq. (I2T!) can be generalized to compute the bias as a function of quasar 
luminosity. Assuming that quasar luminosity closely tracks halo mass, one 
would naively expect an increasing clustering strength for more luminous 
quasars. Whi l e this seems to indeed b e the case for high-redshift quasars 



Francke et al.l (120081 ) ; IShen et al.l (120091 ) , lower redshift surveys fi nd evidence 



instead for a much flatter quasar bias as a function of luminosity. iLidz et al 



( 12006! ) investigated the luminosity dependence of quasar clustering in more 
detail and concluded that by adopting more complicated, non-monotonic 
quasar light curves, such as the ones predicted from numerical simulations 



33 



Hopkins et al.l (j2006al ). a better match to the low-2; data can be found. This 



is because in this picture the Eddington ratio distribution at fixed BH mass 
is much broader, thus mapping both fainter and brighter quasars into similar 
bins of BH/halo mass, and therefore fiattening the bias as a function of 
luminosity. The result of their model is shown in Figure [10] (dotted and solid 
lines), compared to a model (dashed line) which maps quasars into halos 
with a nearly one-to-one relation, and thus predicts a much stronger rise of 
the bias with increasing luminosity. 

Other important constraints can also be derived by studying the scale- 
dependence of the quasar bias. If luminous quasars are predominantly trig- 
gered by merging/interaction events, then their clustering strength at small 
scales should be enhanced with respect to a random non merger popula- 
tion with similar masses and large-scale clustering prop erties. Indeed, such 
a small-scale exce ss has been see n , in b o th low-redsh if t ISerber et al.l (120061 ) 
and high-redshift iHennawi et al.l (|2006f): iMvers et al.l (l2007bl ) quasar popu- 



lations. Interestingly, however, ISerber et al.l (|2006l ) see no such excess in 



the low- redshift Seyfert population, as shown in Figure [TT|, which plots 
the excess small-scale c lustering of quasars and Seyferts with data taken 
from Serber et al.l ( 2006 ) for optical quasars (red filled squares) and Seyferts 
(black open squares). The lines show the predicted predicted bias as a func- 
tion of scale radius determined from the correlation function at each radius 
relati ve to that at large scales, as predicted from t he evolu t ionary mod- 



els of iHopkins and HernquistI (120091 ). Interestingly, iLi et al.l ( 120081 ) using 
SDSS data showed that while specific star formation rates of galaxies are 
higher if they have close neighbors, close neighbors are not associated with 
any similar enhancement of nuclear activity. This leads to conclude that 
star formation induced by a close companion and star formation associated 
with BH accretion are distinct events. They also suggest that these events 
might still be part of the same physical process, for example a merger, pro- 
vi ded they are separat ed in time. These empirical results, as well as those 
by ISerber et al.l (120061 ) , support a scenario in which lower luminosity, often 
lower mass BHs, as usually associated with Seyfert galaxies, might be more 
likely triggered by in-situ phenomena, at variance with the violent dynamical 
event s which are believed to be at the origin of the high- 2:, luminous quasars 



e.g.. ISanders et al.lll988l . see also iKawakatu et al.ll2003l and §|8]). 



Particularly interesting constraints can be derived fro m the recen t ^ > 3 

(l2007h . 



clustering na e asure ments of luminous quasars in SDSS by lShen et al 



White et al.l (120081 ) have applied Eqs. (I2U1) and (1^ to the data, concluding 
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Figure 10: Bias-squared of quasars as a function of their luminosity. The solid and dotted 
lines show the luminosity dependence of quasar clustering in a model with a non-linear 
light curve as predicted from numerical simulations, while the dashed line refers to a 
strictly monoton ic model between quasar luminosity and halo mass. Figure taken from 



Lidz et all (|2006h 
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Figure 11: Excess small-scale clustering of quasars and Seyferts. The data are from 
Serber et al. ( 20061 ) for optical quasars (red filled squares) and Seyferts (black open 



squares). Plotted is the predicted bias as a function of scale radius determined from 
the correlation function at each radius relative to that at large scales (i.e., dividing out 
the best- fit large-scale dark matter correlation function). The dotteddashed line shows a 
constant bias as a function of scale (i.e., no excess of small-scale densities, expected if 
the galaxies are a random subpopulation with no preference for small groups), the solid 
(dashed) orange line shows the predicted excess on small scales expected for recent merger 
populations from cosmological models a nd simulations (reflecting a pr eference for small 
group environments). Figure taken from Hopkins and Hernguist ( 2009() . 
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that the strong clustering measured at z ~ 4 can be understood only if quasar 
duty cyc les are high arid the in trinsic scatter in the luminosity-halo relation 
is small. IShankar et al.l (l2008bl ) took a further step by jointly considering the 
evolution of the BH-halo relation and the BH mass function, as constrained 
by the observed AGN luminosity function and clustering. They examined 
constraints on the host halos, duty cycles, radiative efficiencies, and mean 
Eddington ratios of massive BHsat 2; > 3, as imposed by the clustering 
measurements of IShen et al.l (120071 ) and by a variet y of measurernents of 



the quasar luminosity function at 3 < 2: < 6 (e.g.. iHopkins et al.l. l2007b 



Shankar and Mathuil . 2007; Fontanot et al. . 2007; Shankar et al. . 2009bl . and 



references therein). They linked BH growth to halo growth via a relation 
between BH mass and halo virial velocity V^ir. They found that the strong 
clustering of AGNs observed at z = 3 and, especially, at 2; = 4 implies that 
massive BHs reside in rare, massive dark matter halos. In turn, reproducing 
the observed luminosity function as implied by the growth of dark matter 
halos, requires high efficiency e and/or low Eddington ratio A, with a lower 
limit (based on 2a agreement with the measured z = A correlation length) 
e > 0.7A/(1 + 0.7A), implying e > 0.17 for A > 0.25. The rapid drop in 
the abundance of the massive and rare host halos at 2; > 7 also implies 
a proportionally rapid decline in the number density of luminous quasars, 
much stronger than simple extrapolations of the 2; = 3 — 6 luminosity function 
would predict. 

Figure [T^ shows a compact summary of their main result. The left panel 
shows that only assuming maximal values of the duty cycle A the predicted 
clustering strength (solid and dashed lines) can be, at least marginally, con- 
sistent with the data aX z = 4 (grey bands). The right panel shows with 
horizontal stripes the integrated BH mass density above L = 10"^^ erg s~^ and 
4.0 < z < Q derived from the AGN luminosity function for three different 
values of the radiative efficiency, as labeled. The triple dot-dashed, solid, 
long-dashed, and dot-dashed lines indicate the BH mass density implied by 
the AGN luminosity function at 2; = 4.0 integrated over mass assuming a 
mean Eddington ratio A = 0.1, 0.25, 0.5, 1, respectively, as a function of duty 
cycle A. High duty cycles imply quasars to reside in less numerous, more 
massive and more biased halos. Therefore reproducing the high emissiv- 
ity of luminous quasars requires BHs accreting at a significant fraction of 
the Eddington luminosity (A > 0.25) to radiate at high radiative efficiencies 
(e > 0.15) to match the AGN luminosity function and its evolution with 
redshift. 
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Figure 12: Left panel: predicted z = 4.0 clustering correlation length ro computed above 
the minimum survey sensitivity as a function of d uty cycle, having adop ted the luminos- 



ity function derived from the reference model by IShankar et aL ( 2008b[ ). The solid and 



long-d ashed lines refer to the tq implied by using different bias formulae (see Shankar et al 
2008br ). Dark and hght shaded bands show the la range of the Shen et al. (2007) measure- 
ments at these redshifts. Maximal values of the duty cycle predict a clustering strength 
only marginally consistent with the data at z = 4. Right panel: horizontal stripes show 
the integrated black hole mass density above L = 10'*^ergs~^ and 4.0 < z < 6 derived 
from the AGN luminosity function for three different values of the radiative efficiency, as 
labeled; triple dot- dashed, solid, long-dashed, and dot-dashed lines indicate the black hole 
mass density implied by the AGN luminosity function at z — 4.0 integrated over mass 
assuming a mean Eddington ratio A = 0.1,0.25,0.5, 1, respectively, as a function of duty 
cycle A. High duty cycles require high radiative efficiencies or low Eddington ratios to 
reconcile the cumulative AGN emissivity with the b lack hole mass density in rare halos 
(see text). Figures taken from lShankar et al.l (|2008br ). 



Wyithe and Loebl (120091 ) also suggested that the strong z = 4 clustering 



could be explained by tying qua sar activity to recent l y merged halos, which 

migh t have stronger bias (see also Wechsler et al.|[2006 ; Furlanetto and Kamionkowski 



20061 : IWetzel et al.l l2009l ). However, halos with substantial "excess" bias 



might be too rare to satisfy duty cycle constraints, and/or the amount of ex- 
tra bias mig ht be too modest to ma ke any significant difference in the models 
described by lShankar et al.l (j2008bl ). Only numerical simulations will be able 
to provide significant answers in this direction. 
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8. AN OVERVIEW ON COEVOLUTION MODELS OF BLACK 
HOLES AND GALAXIES 



?. i. GALAXY DICHOTOMY AND AGN FEEDBAGK 

Galaxies can be classified into two families, which reflect distinct morpho- 



l ogica l, chemical and evo l ution ary fea tures. Following, e.g., iKauffmann et al 



( 120031 ). iHeckman et al.l (120041 ). and iDekel and Birnboinj (120061 ). here are 
summarized the main pieces of observational evidence supporting this bi- 
modal classification. There is a gap at a color of u — r ~ 2, where galaxies 
naturally divide into a blue and a red sequence. The former are on av- 
erage younger, less massive and disky {late-type galaxies), while the latter 
are more massive older and boxy, spheroidal {early-type galaxies). Blue, late- 
type galaxies dominate the mass function below the knee, L*, of the Schechter 
function while red, early-type galaxies take over above L*. 

There is a remarkable correlation between the colors and luminosities of 
elliptical galaxies, with brighter galaxies being redder. Values of the ele- 
mental ratios [a/F e], in particular for Mg, have been found to incre ase with 
galactic mass (e.g., Terlevich and Forb ei. l2002l : iBernardi etaP . l2003h , imply- 
ing that the present rate of Supernovae (SN) I is ten times higher than the 
rate of SNII. The latter is close to being a tracer of the instantaneous star 
formation rate. These observations thereby constitute significant evidence 
that massive galaxies formed most of their stars at 2; > 2-3 in a relatively 
short burst. Direct evidence that massive galaxies with Mgtar > 10^^ Mq were 
in pl ace at z > 2, is provided by recent -ft"-ba nd and submillimeter surveys 

(e.g., 
lation 



Cimatti et all, I2OO2I: ICreve et all IgOOsh . The fundamental plane re- 



2009 



(e.g. 



Djorgovski and Davisl . 1987 : Treu . 2004; Shankar and Bernardi . 
and references therein) observed up to 2; ~ 1, is also consistent with 
massive spheroids being old and quiescent. 

According to the standard cosmological paradigm of structure forma- 
tion and evolution, dark matter halos have grown hierarchically, through the 
merging together of smaller units into ever larger systems. In this scenario, 
galaxies form inside this hierarchically growing system of halos. However, 
hierarchical models do n ot predict the correct abundances of galaxies at high 



and low redshifts (e.g., iFontana et al.l . l2006l : iMarchesini and van Dokkunj . 



20071 ) and, at variance with observations, they also tend to generate negative 



gradients for the [a / Fe] abundance with respect to stellar mass, as more Fe 
is formed in more massive systems whi ch comp l ete th eir formation at later 



times in pure hierarchical models (e.g.. iThomad . Il999l ) 
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Therefore in recent years theoretical models of galaxy formation have in- 
serted AGN feedback as a fundamental ingredient for successfully matching 

several pieces of i ndependent data, from the stellar mass function to chemical 

abunda nces (e.g.JWvithe arid Loebl.l2003l:lGranato et al.l. 12004: IScannapieco and Oh 



2004: 



Pi Matteo et al.l 



2005 



Sazonov et al 



2005: Vittorini et al 



2005 



Bower et al 



20061: Cattaneo et al.l.l2006l:ICroton et al.l.l2006l:lDe Lucia et al.l.l2006l:lHopkins et al. . 



2006al:lLaDi et al.l.l2006l:IShankar et al.l.l2006l:lMalbon et al.l.l2007[ 



20081 : iBenson et all . l2007l : iMonaco et aP . l2007l : ICook et all . l2009h . How 



Menci et al 



ever, clear and statistically significant evidence for AG N feedback in galax- 



ies is still lacking, despite some significant studies fe.g.. IChartas et al.l.l2003l: 



2008cl . 



a 



King and Pounds! . l2003l : ISchawinski et al.l . 120071 : iDai et al.l . l2008l : IShankar et al. 



AGN feedback consists of injection of energy and momentum acting on 
the gas and dust of the host galaxy. Its implementation in theoretical mod- 
els differs significantly between different research groups; the energy can be 
in fact injected into the system in a radiative, thermal, or kinetic mode, 
and/or at different times during the evolution of the system, and/or under 
different conditions. These models are, nevertheless, similar in the aims of 
removing gas from the galaxy, quenching star formation, and turning the 
galaxy from blue to "red and dead". AGN feedback has also been proven 
to predict the correct slopes and normalizations of th e observed sca l ing re - 



lations between BHs mass and galaxy properties (e.g., I Silk and Reed . Il998l ). 
such as the stellar dispersion velocity a or ii'-band luminosity, although feed- 
back is not a necessary condition to s ubstantiate such scaling relations (e.g.. 



Miralda-Escude and KoUmeierl . l2005l ). 



The physics of accretion onto the central BH remains, however, poorly 
understood. The problem of how to remove angular momentum from gas 
in the inner regions of galaxies allowi ng it to accrete o n to the central BH, 
is still far from settled. Furthermore, [Thompson et al.l (120051 ) pointed out 
that even if gas is efficiently transferred to the central regions through some 
global torque, then the high densities will lead all gas to be easily converted 
into stars at the 1-10 pc scale, leaving no fuel for the BH. 

Most SAMs now include some basic recipes to account for BH growth 
and the co-evolution with their host galaxies, although the assumptions and 
results sometimes differ significantly among different groups. Some models 
predict, for example, a merger dominated growth of BHs at la te times, while 
others state that gas accretion always dominates the evolution. iMalbon et al. 
(120071 ). within the context of their hierarchical model of galaxy evolution. 
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find that while lower mass BHs mainly grow through gas accretion during 
starbursts, more massive BHs grow first by accretion at higher redshifts, and 
then predominantly by merging of pre-existing BHs at lower redshifts. This 
conclusion is not shared by other SAMs, wh ich instead grow all BHs mainly 
through gas accretion (e.g. 



Lapi et al.l 120061 ). 



SAMs have adopted several different prescriptions for triggering accre- 



tion onto the central BHs, such as radiatio n -drag (e.g., iKawakatu et al. 



20031 ) , bar-instabilities fe.g.. iBegelman et al.l . Il984j ). cloud-cloud collisions 



e.g.. iHopkins and Hernguistl. 120061). merging or galaxy interact ions (e.g.^ 

Cavaliere and Vittorini . 20001 : Di Matteo et al. . 2005 : Croton et al. ^2006: H opkins and Hernquist . 



20061 : iHopkins et al.l. l2006al). surface density instabilities fe.g.. ICenl. l2007l). 



stella r feeding (e.g.. lMiralda-Escude and Kollmeieii . l2005l : ICiotti and Ostriker 
20071 ). Nevertheless, most SAMs provide reasonable matches to the same ob- 
servables, such as the optical quasar luminosity function, although sometimes 
their input BH physics might differ significantly. 

While independent, fundamental tests for BH evolution will probably 
need to wait for future gravitational wave experiments, such as LISA (the 
"Laser Interferometer Space Antenna", a joint ESA/NASA mission), impor- 
tant constraints can be set on BH evolution in ways independent of specific 
models, via basic accretion and merging models which also include cluster- 
ing, as already extensively discussed in the previous sections. In this section 
we instead briefiy review some of the main results achieved by several SAMs 
which attempted to reproduce the overall BH evolution. 

8.2. UNDERSTANDING THE EVOLUTION AND SHAPE OF THE AGN 
L UMINOSITY FUNGTION 

Th e two power-law observed AGN bolometric LF (e.g., "i lShankar et al. 



2009bl . and references therein) is an essential tool to probe BH evolution. In 
the overall, a significant decline in the AGN luminous density is observed, 
with the faint-end rapidly steepening at redshifts z < 1.2. 

Such an evolution could be described by two extreme scenarios. On the 
one hand, less massive BHs form continuously over time, to track the popula- 
tion decline to lower luminosities while retai ning nearly constant (o r slightly 
declining) Eddington ratios A = L/M,, as in lKoUmeier et al.l (120061 ). On the 
other hand, the strong AGN luminosity density redshift evolution could also 
be explained by a distribution of already massive BHs shining with ever lower 
luminosities over time, inevitably causing a significant decrease in the overall 
Eddington ratio distribution. From the simultaneous match to the quasar 
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clustering and luminosity function, IWyithe and Loebl (120051 ) concluded that 
bright quasars at z < 2 have a lower abundance than expected from the 
number of merging halos, and discussed how a viable solution would be for 
quasars to shine at sub-Eddington luminosities. Similar questions arise re- 
garding the steepening of the faint-end of the AGN LF: is it produced by 
lower-mass BHs shining at the Eddington limit or is it mostly due to high- 
mass BHs radiating at sub-Eddington luminosities? 

Within the frame work of their numerically-predicted BH light curves, 
Hopkins et al.l (j2006bj) found that the bright-end slope and break luminosity 
time behavior in the AGN luminosity function, directly reflects the formation 
rates of BHs shining close to their peak emission following a major merger, 
while the main contribution to the faint-end is from BHs shining at luminosi- 
ties far below the peak. This scenario would imply a continuous formation 
of massive BHs down to very low redshifts driven by major mergers. On 
other grounds, however, extending major mergers as main triggers for the 
overa ll quasar popula t ion at all times, might g i ve ris e to several problems 



e.g., iBournaud et al.l . 120071 : iNaab and Ostrikerl . l2009l ). 



8.3. THE HIGH-REDSHIFT EVOLUTION: FAST GROWTH OF MAS- 
SIVE GALAXIES AND BLACK HOLES 

At high-2;, the main SAMs fail to grow BHs fast and massive enough 
to reproduce the relatively high n umber of lu r ninous quasars observed at 
high- 2; and their extreme clustering IShen et al.l (120071 ). High-resolution nu- 
meric al simulatio ns have shown that DM halos grow in a two-phase mode 
(e.g., Zhao et al. 1, [2003). The inner core and central potential of the host 
halo is built during a fast, chaotic, merging phase of multiple subhalos 
after which a longer, smoother accretion phase occurs. During the lat- 
ter phase, mass is preferentially added to the outskirts of the halo with- 
out affecting the inner regions much. This behavior may help in preserving 
the scaling relations between the central BH and host galaxy, and, in gen- 
eral, hmit the role of environment in th e formation of luminous galaxies 



e.g., iBournaud et all . l2007l : iTinker et all . l2008l : ICook et all . 120091 ). In the 



DM potential wells, SAMs attempt to speed up the formation of massive 
galaxies and their BHs by allowing baryons to cool in high clumping fac- 



tors 



e.g. , iGranato et al.l . |2004J : iMaller and BuUockl . |2004J : IScannapieco et al 



Dekel and Birnboim 



2005 ) and/or cold streams flowing directly to the center (e.g. 
20061 ). For example, by using enhanced clumping factors and short forma- 



tion timescales, Granato et al. (2004), ICirasuolo et al.l (120051 ) . Lapi et al 
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(2006) and Shankar et al. (2006), were able to reproduce several dynami- 
cal, photometric, and statistical properties of BHs, spheroidal galaxies and 
AGNs. 

At the same time, SAMs need to cope with the difficulty of reproducing 
the non-linear metallicity behaviors observed in early-type galaxies, with the 
more massive galaxies being [a /Fe] -enhanced with respect to the less massive 
ones. The Fe production is roughly proportional to the (delayed) Supernova 
la rate. Thus, reducing the burst timescale in the more massive galaxies, 
prop ortionally reduce s their amount of Fe inducing an [a/Fe]-enhancement 
(e.g.. lMatteuccilll994j ). In practice, reproducing, and physically explaining, 
this differential trend in metallicity is a challenge f or SAMs. Gr a nato et al. 
(2004), similarly to o ther models s uch a s those of ICroton et al.l (120061 ) . De 



Lucia et al. (2006), iPipino et al.l (120091 ) . naturally find the more massive 
galaxies to be, on average, the oldest and most metal rich. 

However, to get the right [a /Fe] -enhancement some extreme assumptions 
need to be made. Granato et al. (2004), for example, allowed for Super- 
Eddington accretion at high redshifts ( this assumption has a lso been used 



in the recent numerical simulations by iDi Matteo et al.ll2008l ). thus boost- 



ing BH growth which more rapidly meet the condition for self-regulation 
and quench star formation earlier than less massive galaxies. iGranato et al. 



( 120061 ) successfully adop ted the same prescrip tions to reproduce the bulk of 
the 850/i counts (see also lFontanot et al.ll2009l ) observed with the submillime- 
ter a rray SCUBA, a difficult task for hierarchical galaxy formation models 



e.g., iGranato et al.l . l2000l ). By allowing for substa ntial BH growth duriii g 



the dust-obscured phase triggered by radiation drag (IKawakatu et al.ll2003l ). 
they also predicted the correct number co unts of AGN sources era itting above 
a certain X-ray flux limit, as measured by [Alexander et al.l (120051 ). Lapi et al. 
(2006) using the same model to grow BHs to at z > 2, and a significantly large 
scatter in the quasar luminosity-halo relation, were then able to reproduce 
the bright end of the observed AGN luminosity function (see the previous 
Authors for details). However, a large scatter in the luminosity-halo relation 
of high-z quasars might be difficult to reconcile with the extreme clustering 
signal measured for those quasars (see § [7]), and more secure data and more 
detailed modeling is required to set final conclusions. 

Similar conclusions on the growth of the most massive BHs at high-^ 
were reached by MaruUi et al. (2008; similar results were found by Fontanot 
et al. 2006 using the MORGANA model). Although their semi-analytical 
framework, developed on the outputs of the Millennium Simulation, better 
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tracks, in principle, DM merger histories than the analytical treatments by 
Lapi et al. (2006), their final conclusions on the baryonic physics is similar. In 
fact, while they can reproduce the observed AGN luminosity function at low 
and intermediate redshifts, aX z > 1 they under-predict the number density 
of bright AGN, regardless of the BH mass accretion rate and light curve 
model assumed for each quasar episode. Despite the considerable freedom in 
choosing the Eddington ratio distributions, they failed to find a model able 
to simultaneously match the observed BH scaling relations, the BH mass 
function, and the AGN luminosity function, especially at high redshifts. Also 
allowing for larger seed BH masses, they were still not able to fit the high-^ 
luminosity function. They interpreted this failure as an indication that the 
present theoretical framework is itself inadequate to account fo r the full AGN 



phenomenon. Analogously to Lapi et al. and Granato et al., iMaruUi et al 



( 120081 ) also suggested that a significant improvement in the model might arise 
by assuming an accretion efficiency that increases with redshift, for example 
by forming more con centrated disks in the past, making more efficient the 



BH feeding (see, e.g.. lMo et aPllQQSh 



8. 4 . THE INTERMEDIA TE REDSHIFT E VOL UTION: INTER A CTIONS, 
BLACK HOLES RE- ACTIVATIONS, ACTIVITY IN DISKS 
Although efficient, the fast accretion phase occurring during the epoch of 
rapid merger events in the high-redshift Universe cannot be the only viable 
formation mode to fully grow all local BHs. For z < 1.5 major mergers 
become much rarer at galactic scales; fewer new massive BHs are then ex- 
pected to be formed. The prevailing dynamical events that could still trigger 
star formation in a galaxy and possible gas accretion onto the central BH, are 
best described as interactions between galaxies, occurring in the small, dense 
groups that begin to form at these epochs. The accretion is also expected to 
be supply-limited, due to progressive dilution of the host gas reservoir and 
lack of fresh gas resuppl y from gas-rich mergers. 



As discussed in (e.g., ICavaliere and Vittorinil . l2000l : iHopkins et al.l . l2008l . 



and references therein), small galactic groups with virial mass > 10 Mq 
could provide particularly suitable sites for the hosts to interact with their 
companions. An interaction of the host with a group companion will perturb 
the galactic potential, and destabilize a fraction / of the cold gas mass m{z) 
in the host from kpc scales to the inner regions, possibly tri ggering a cir cum- 



nuclear starburst (see, e.g., the numerical sina ulations by iMihod Il999l . and 



the IR observations by iFranceschini et al.lll993l ) and reactivating the central 
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BH to a level of luminosity L oc fm{z). Direct evidences of star-formation 
activity at intermediate redsliifts connected with c learly interact ir ig gal axies 



have now been provided by several groups (e.g., lUrrutia et al.L l2008l . and 
references therein), who also find that a large fraction of AGNs have close 
companions or show signs of ongoing or recent interactions, although locally 
and/or at lower luminosities this association between AGN activity and close 
co mpanions might break down, as discussed in § [71 



Vittorini et al.l (120051 ) have carried out a detailed semi-analytical com- 



putation of interactions and minor /major merging of galaxies within the 
framework of the hierarchical growth of dark matter halos. They found that 
the strong luminosity decline with redshift of the AGN luminosity function, 
given by L oc / X m{z), is fully accounted for by the continuous exhaustion 
of gas m{z) driven by galaxy interactions. The highly stochastic nature of 
the accretion events described by interactions, parameterized in the factor 
/, predicts a largely increasing scatter in the Eddington ratio distributions 
at lower redshifts, peaked at lower and lower accretion rates, in line with the 
several recent empirical findings discussed in § [3 

Other possible ways to reproduce the bright end o f the quasar luminos- 



i ty fu nction have been extensi vely discussed in, e.g., IScannapieco and Oh 
( 120041 ). and lCroton et al.l (120061 ). These models mimic the fall-off and steep- 
ening of the quasar luminosity function at low redshifts by assuming that 
late AGN feedback in massive halos prevents cooling fl ows in massive halos 
(the s o-called "radio-mode" feedback). More recently, ICavaliere and Menci 
( 120071 ) have presented the full predictions of their semi-analytic model for 
the co-evolution of BHs and galaxies, which includes galaxy interactions, 
among other dynamical triggers. Their model predicts a "bimodal" Edding- 
ton ratio distribution, with a peak around A ~ 0.01, mainly contributed by 
low-luminosity AGNs in red galaxies, and a peak at A < 1 contributed by 
lower mass BHs in bluer, star-forming galaxies. 

It is c lear, especially after the significant results obtained from SPSS 



data (e.g., Kauffmann et al.l . 120031 : iKauffmann and Heckmanl . 120081 : iLi et al. 



20061 . I2OO8I ). that a significant part of the AGN activity in the local Uni- 



verse must be triggered by other mechanisms other than major mergers, 
such as bar instabilities, radiation drag, stellar winds. Moreover, observa- 
tions have un yeiled the presence of AGNs also in isolated disk-dominated 



galaxies (e.g., ISatyapal et al.l . 120081 ). as discussed in § 13. 1[ AGN feedback 



migh t have been possibly efficient du e to the lower mass in the central BHs 



e.g., iKawakatu and Umemural . l2004 IShankar et al.l . 120061 ). and the metal 
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rich gas outflow may be expected to recoUapse towards the center on a rela- 
tively short timescale, due to the shallo wer potenti a l well s that characterize 



the halos hosting disk galaxies today. iTsujimotd (120071 ) finds that if one 
includes a double infall of primordial and enriched gas in chemical evolu- 
tion models of the Galaxy, with the enriched outfiow occurring after 1 Gyr 
since the onset of the star formation, it is possible to explain the observed 
a-element abundances and their non-linear behavior in the Galaxy as a func- 
tion of iron abundance. This model can also simultaneously reproduce the 
sharply peaked metal abundance of disk stars against iron abundance. 

To sum up, most models agree in associating the luminous, high- 2; quasars 
with fast, gas-rich merging events, which rapidly built up the host potential 
wells and the scaling relations discussed in § HJ A longer and milder activity 
phase has probably occurred at lower redshifts, during which in-situ and 
milder dynamical events have regulated the late evolution in the AGN/BH 
population. However, as discussed in this section, many details of this basic 
picture are still missing, and are currently the subject of intense investigation. 



9. Conclusions 

This Review addresses the main topics and aims concerning the evolution 
of the BH population, and discusses several ways to constrain it in ways 
independent of specific models of galaxy evolution. The main results can be 
listed as follows. 

• The local BHMF can be now rather accurately assessed. In more detail, 
the BHMF hosted in early-type galaxies can be obtained on exploiting 
the velocity dispersion or luminosity functions of host galaxies, coupled 
with the Mbu-ct or M-BR-Lsph relationships, respectively, yielding sim- 
ilar results within ~ 30% accuracy. The contribution from BHs hosted 
by late-type galaxies is more uncertain, and is mostly confined to the 
low mass end of the BHMF. The overall BH mass density amounts to 
PBH = (3 - 5) X lO^Mo/Mpc^^ (§ ED. 

• The match both between the local BH mass density and the one pre- 
dicted by integrating the bolometri AGN luminosity function (Soltan 
1982) yields interesting constraints for a mass to radiation conversion 
efficiency e < 0.1, although systematic uncertainties on bolometric cor- 
rections, luminosity functions, and the local BHMF prevent more pre- 
cise estimates (§11]). 
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• The average accretion history of BHs tracks the average SFR, with a 
proportionahty constant of the order of the local one. Several works 
have empirically proven that the ratio between the degree of evolution 
in the local relations is relatively mild, suggesting that the bulk of the 
BH and host galaxy populations should have mostly co-evolved, at least 
up to 2; < 2 (§11]). 

• The shape of the local BHMF yields additional interesting constraints 
on the typical Eddington ratio A ~ 0.5. The alternative possibility 
that most of the mass locked in BHs in the local Universe has been 
accumulated by "dark" accretion (i.e., accretion undetectable by either 
optical or hard X-ray surveys, as in the case of BH coalescence), is 
severely constrained by the above results. In order to make room for 
this possibility one has to assume that the radiative efficiency during 
the visible AGN phases is at the theoretical maximum of e ~ 0.3- 
0.4. But even in this case (unless the bolometric correction is far lower 
than currently estimated) the contribution of radiative accretion to the 
local BH mass density is > 25%, and one is left with the problem of 
fine tuning the radiative and non-radiative contributions in order not 
to break down the match with the local BHMF obtained with radiative 
accretion alone. One would also face the problem of accounting for the 
tight relationships between BH mass and mass or velocity dispersion 
of the spheroidal host, naturally explained by feedback associated to 
radiative accretion. The analysis of the accretion history highlights that 
the most massive BHs (associated to bright optical quasars) accreted 
their mass faster and at higher redshifts (typically at z > 1.5), while 
the lower mass BHs have mostly grown at z < 1.5 (§§ [5] and [6]). 

• Quasar clustering provides additional, independent constraints on the 
main parameters governing the BH evolution. It has been shown how 
AGN clustering can set important constraint on the possible redshift 
and/or mass dependence of the typical A and e, thus helping us in dis- 
criminating among a number of successful models. Clustering is in fact 
a direct measure of the masses, and therefore number densities, of the 
halos hosting the quasars. Matching the quasar host number density 
with the one derived from accretion models (Eq. [H]), in turn yields 
interesting constraints on radiative efficiency and Eddington ratios as 
a function of time and/or BH mass (§ [7]). 
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• Most galaxy formation models now include some basic recipes to ac- 
count for BH growth and the coevolution with their host galaxies, al- 
though it has been discussed in this Review that the assumptions and 
results sometimes differ significantly among different groups (§[8]). The 
model-independent results discussed above will help shedding light on 
the true physical mechanisms which regulate BH growth across time. 
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